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Abstract. A modified version of the double potential formalism for the elec- 
trodynamics of dyons is constructed. Besides the two vector potentials, this man- 
ifestly duality invariant formulation involves four additional potentials, scalar 
potentials which appear as Lagrange multipliers for the electric and magnetic 
Gauss constraints and potentials for the longitudinal electric and magnetic fields. 
In this framework, a static dyon appears as a Coulomb-like solution without 
string singularities. Dirac strings are needed only for the Lorentz force law, not 
for Maxwell's equations. The magnetic charge no longer appears as a topologi- 
cal conservation law but as a surface integral on a par with electric charge. The 
theory is generalized to curved space. As in flat space, the string singularities 
of dyonic black holes are resolved. As a consequence all singularities are pro- 
tected by the horizon and the thermodynamics is shown to follow from standard 
arguments in the grand canonical ensemble. 
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1 Introduction 

Reissner-Nordstr0m black holes with both electric and magnetic charge 

ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 (d6 2 + sin 2 6d<p 2 ) } 

I 2M Q 2 + P 2 
N=\l + V 2 , (1.1) 

A = -9.dt + P(l -cos9)d0, (1.2) 
r 

are generally excluded in a discussion of uniqueness theorems and geometric derivations 
of the first law because the gauge potential is singular along a string that intersects the 
horizon and goes to infinity 03 |2j. Exceptions can be found in for stationary and 
axisymmetric perturbations and in flU where dipole charge contributions to the first law 
for five dimensional black ring solutions are investigated by dealing directly with diver- 
gent potentials on the horizon. Dyonic solutions were also excluded in an investigation 
of duality of electric and magnetic black holes using Euclidean methods [|5]|. 

Nevertheless, for variations of the three parameters the first law 



6M = -£-6A + <f> B SQ + il>H6P, (1.3) 



•57T 



where 



A = M 2 — (Q 2 + P 2 ), r± = M±VA (1.4) 
' A = 8n[M 2 - Q2 + p2 +MVA], (1.5) 



Kj . , . ........ 

2r\ ' 2 
Q P 

<pH = — , Wh = —, (1-6) 

can easily be inferred from the purely electric case by using a duality argument. Further- 
more electric-magnetic black hole duality has been extended to the case of dyons in the 
canonical ensemble by using the manifestly duality invariant double potential formalism 
[0. In its original version 0, this formalism involves as dynamical degrees of freedom 
two vector potentials. An independent rederivation [8 J has been written with two addi- 
tional scalar potentials which are spurious because they appear only as a part of a total 
derivative of the action. In the black hole context O, coupling to external static sources 
can be made either through fixed strings, spherically symmetric nondynamical longitudi- 
nal fields, or intermediate combinations. Finally, the coupling to dynamical dyons with 
the help of dynamical strings has been studied in detail in 0EQ1, including a proof of 
equivalence with Dirac's original theory lfTTl[T2ll and a derivation of the appropriate quan- 
tization condition [fT3l[T4ll . 

What we will do in this paper is introduce potentials for longitudinal components of 
electric and magnetic fields. This has the effect of making the two scalar potentials non- 
spurious as they now appear as the Lagrange multipliers for the divergence constraints on 
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electric and magnetic fields. We thus increase the redundancy of the description in such 
a way as to have twice as much gauge invariance as in standard Maxwell theory. 

Now, taking into account all the results described above, this extension of the double 
potential formalism is rather straightforward and seems hardly worth the effort. We beg 
to differ. 

First of all, the electric and magnetic potentials produced by a static dyon both appear 
as Coulomb-like solutions in a single, manifestly duality invariant formulation without 
any stringlike singularities. In this framework Dirac strings are only needed in order 
to produce the correct Lorentz force law from an action principle for dynamical point- 
particle dyons. 

In curved space, the new formulation is ideally suited for the description of black hole 
dyons. As in flat space, their string singularityis resolved and a geometric derivation of 
the first law can be done along standard lines because all singularities are now protected 
by the horizon. This is a direct consequence of the intriguing transmutation into a surface 
integral of the magnetic charge which appears as a topological conservation law in the 
standard approach. Since there is no quantization condition on magnetic or on electric 
charge for a single dyon and because of the presence in the formalism of both chemical 
potentials, thermodynamics and Euclidean computations can be performed in the grand 
canonical ensemble, thus circumventing arguments of fl2l[5]|. 

In the next section, we discuss our formulation in Minkowski space in the case of 
fixed external sources. Section|3]is devoted to Dirac strings and dynamical point-particle 
dyons. We finally write down and analyze the appropriate action for curved space and 
discuss applications in the context of black hole physics. 

2 Extended double potential formalism in flat space 

In this section we present an action principle for electromagnetism in the presence of elec- 
tric and magnetic sources which is manifestly duality invariant. Both electric and mag- 
netic Gauss constraints are dynamical and appear in the action with their corresponding 
Lagrange multipliers. For a static dyon, the solution of the field equations is Coulomb- 
like, both in the electric and the magnetic sector. We show that the theory can be gauge 
fixed so as to coincide with standard electromagnetism and conclude the section by show- 
ing that Lorentz invariance, while not manifest, is nevertheless realized through canonical 
generators very much as in the standard Hamiltonian formulation of electromagnetism. 
This suggests, as we will explicitly show in the last section, that the theory can be gener- 
alized to curved space. 
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2.1 Action, duality, and gauge symmetries 

The dynamical fields of the theory are A^, C a , a = 1, 2. Here, A^ = (A^, Z M ) are the 
standard and new potentials. The additional fields C a = (C, Y) make up the longitudinal 
parts of magnetic and electric fields B a = (B, E) according to 

B a = V x A a + VC a . (2.1) 

The external magnetic and electric currents j aM = (k^,j^) are conserved, d^j^ = 0. In 
this section, we assume that they correspond to the currents produced by a single point- 
particle dyon. We consider the action 

I[A^C a ] = I M [Al,C a ]+Ij[A^j a % (2.2) 

where 

I M [A«, C a } = \ j d 4 x [e ab (B a + VC a ) ■ {d^ - VA b ) - B a ■ B a j , (2.3) 

is the substitute for the usual Maxwell action and 

h[A%r} = J d 4 xe ab A a ^ (2.4) 

is the "interaction" action. Here e ab is skew-symmetric with e 12 = 1, and indices a,b, . . . 
raised and lowered with the Kronecker delta. The action (12.21) is manifestly invariant 
under simultaneous duality rotations on (A^, C a ; j aM ) 

5 D A° = e ab A b „ 6 D C a = e ab C b , 5 D f» = e ab j£. (2.5) 

It is also gauge invariant under 

6 X A« = d,X a , 5 x C a = 0. (2.6) 



2.2 Equations of motion and point-particle dyon 

The Euler-Lagrange equations of motion associated with (|2.3I) are easily shown to be 
equivalent to Maxwell's equation with magnetic and electric currents. Indeed, variations 
with respect to Aq give the constraints 

V • B a = S7 2 C a = j 0a . (2.7) 

Variations with respect to C a imply the equations 

V 2 C a = e ab (V-d A b -V 2 A b ). (2.8) 
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The fields A®, C a are auxiliary in the sense that, under suitably boundary conditions at 
spatial infinity, their equations of motion can be solved for A®, C a in terms of all other 
fields, without the need for intial conditions. 

Variations with respect to A a yield Maxwell's equations in the form 

- e ab d B b + VxB a = t ab f. (2.9) 



As a consequence, if the electromagnetic field tensor F is expressed in the usual way 
in terms of electric and magnetic fields, F 0i = —Bf, = 6ijkB lk , it follows that both 
dF and d*F vanish outside of sources on account of the Euler-Lagrange equations of 
motion. 

In the case of a single point-particle dyon at the origin with charges Q a = (P, Q), for 
example, 

j a »{x)=A*Q a 8%8\x), (2.10) 
instead of (11.21) . Maxwell's equations in the above form are now solved by 

A « = -^Ldt, C a = -^-. (2.11) 

This solution resolves the string-singularity of the standard formulation. It is unique in 
the transverse gauge V ■ A a = with vanishing boundary conditions on A a , C a . 



2.3 Canonical structure and degrees of freedom 

By using integrations by parts and decomposing A a = A aT +VM a , with M a = (Ma,Mz), 
the free action (12.31 ) can be written in the form 



I M [A aT ,A a ,M a ,C a ] = J d A x 



V x Z T ■ d A T + V 2 Yd M A - V 2 Cd M z - 

(2.12) 



~\f ■ E — \b ■ B — A V 2 Y + Z V 2 C 



2 2 



where E = V x Z + VF, B = V x A + VC. This shows that the canonically conjugate 
pairs are (A T , -V x Z T ), (V 2 Y, M A ) and (-V 2 C, M z ) so that there are 4 conjugate 
pairs per spacetime point. 

Variation with respect to the Lagrange multipliers Z imposes the first class constraint 
V 2 C = 0. Partial gauge fixing to the standard co variant description can be achieved 
by requiring the longitudinal part of the second vector potentials to vanish, M z = 0, 
and gives back the usual Hamiltonian description of electromagnetism. Complete gauge 
fixation is then achieved, as usual, by solving the electric Gauss constraint V 2 F = 
associated with the Lagrange multiplier A together with the gauge condition Ma = 
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0. The gauge fixed theory contains 2 physical degrees of freedom per spacetime point 
described by the transverse vector potential A T and its canonically conjugate variable 
-E T = -V x Z T , as it should. 

For later use, we note that 

{A m {x), B b i(x')} = -e ab 5 lj 5 3 {x, x'), {C a {x), B bj (x')} = 0, 
{M a (x),C b (x')} = e ab \7- 2 5 3 (x,x'), {B m (x), B bj (x')} = e ab e ijk d k 5 3 (x, x'). (2.13) 



2.4 Duality, gauge, and Poincare generators 

The Hamiltonian and constraints associated with the first order action I M [A a , C a ] are 

H= [ d 3 x I B a ■ B\ g a = e ab V ■ B b , (2.14) 



The duality generator is the 5*0(2) Chern-Simons term [7] suitably extended to the lon- 
gitudinal potentials, 

D = -\ J d 3 x (B a + VC a ) ■ A a . (2.15) 

It commutes with the Hamiltonian and the other Poincare generators introduced below, 
but is only weakly gauge invariant, 

{g aj D} = e ab g b . (2.16) 

The duality transformations (12.51 ) on the canonical variables Af, C a are generated through 
SdAI = {A®, D}, 5oC a = {C a , D}. The extension to the Lagrange multipliers is dic- 
tated by (12.151) and the requirement that the first order action I M [A\ C a ] is invariant. In 
the same way, the gauge transformations 5\ in (12.61) are generated by 

T[A] = J d 3 xg a \ a . (2.17) 

In this expression, the generators are smeared with the arbitrary functions A a defining the 
gauge transformation in (12.61) . 

A general Poincare generator may be written as 

T(u,a) = \uo^ -a„P",. (2.18) 

Here J^ v and are the individual Poincare generators, and cu^, a M the corresponding 
parameters defining the transformation. The generator of time translations is the Hamil- 
tonian, P° = H. The Lorentz generators may be decomposed as 

\u^ v = \uo^ k J k - oo i0 K\ (2.19) 
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The Poincare generators are related to the symmetric energy-momentum tensor with com- 
plete electric and magnetic fields as follows: 



T oo = 1 ^2 + T io = (g x gy t (2.20) 

P> = f*,T*, /- = -/A(^-^), (2.2,, 



where J' k = J k and K' — K — x°P 



d 3 x e nh B a xB\ J = I d 3 x e nh B a (x ■ B b ) 



K = J d 3 xx(^B a ■ B a ). (2.22) 

One can then show by direct computation of the Poisson brackets that these generators 
form a representation of the Poincare algebra, up to terms involving the constraints. We 
will prove this explicitly in Sec. 14.31 and show that 



where 



{T(u, a),T(0, &)} = T([u, 9},ub- da) + T[[f , V } B ], (2.23) 



77(0,6)" = -O^^ + V 1 , 
[Z,r)} a B = B^e^rf-e^B^e^-rf?). (2.24) 



If we then define 



a[ = di + uj 0i x° (2.25) 



and all other parameters unchanged, the conserved Noether charges generating the Poincare 
transformations as canonical transformations on the fields are 

Q(u,a) =T(u,a'). (2.26) 

Indeed, deriving (12.231) in terms of a' (a) with respect to b and putting 9 = gives 

{H,Q(u,a)} = §- t Q(u,a) - J d 3 x g a e ab B u C(u,a'). (2.27) 

As a consequence, the Noether charges are conserved on the constraint surface, as they 
should, and the Poincare transformations of the canonical variables, 5qC°" = {C a , Q} 
and 5qA°; = {A%, Q}, can be extended to the Lagrange multipliers so as to leave the 
action invariant. Explicitly, with the understanding that £ = £(u, a' (a)), 

5 Q C a = 0, (2.28) 

5 Q A« = d^Q-e^B^-e^B^, (2.29) 

6 Q B? = -e« k d j ((? b B bk ?)-d j (B a SC)+dj(B ai ?), (2.30) 

5 Q A a = d \ a Q + e ab B bl C, (2.31) 
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where 

A« = -e a6 V- 2 9 i (S*e°) + V-%(e^ k B^ k ). (2.32) 



3 Dynamical point-particle dyons 

We show in this section that for sources that correspond to dynamical point-particle dyons, 
a consistent action principle that makes the dyons evolve according to the Lorentz force 
law needs Dirac-type strings and requires a veto giving rise to the standard quantization 
condition. We then show equivalence with Dirac's original, manifestly Lorentz invariant 
formulation. 



3.1 Dyons and Dirac strings 

We begin by reviewing the use of Dirac strings in the theory of magnetic monopoles. 
Let us first fix conventions. Define e ai ...a„ = e ai "" a ™ to be totally skew- symmetric with 
ei... n = 1. The Levi-Civita tensor is £ ai ...a n = \Ag~\ e ai...a n - Indices on this tensor are 



raised with the metric, which implies that e ai - a ™ = i^e" 1 -"" where a is the signature 

vial 



of the metric. Our convention for the dual is (*uj p ) ai ...a n ^ = ^ bl '" bp e bl ... bp a 1 ...a n ^ p - 

Consider a (d+l)-dimensional surface in flat 4-dimensional spacetime parame- 
terized by (r, ai,...,a d ), 

x » = v i*(t, a 1 ,...,a d ). 
Associated with this surface, define the d + 1 form H^ d+1 with contravariant components 

H^^ d+1 (x) = J 5 {A \x-v)dv^ A ••• Adv^+K (3.1) 
It is straightforward to show that if <9£ is the boundary of E, then, 

d*#s d+1 = *H dT , d+1 . (3.2) 

In the Dirac theory, the worldline T : x^ = z^(r) of a magnetic pole of charge g defines 
the magnetic current 

JL 9 = 9H£- (3-3) 

The worldline is the boundary of the worldsheet of a Dirac string £ : x M = y M (r, a). 
Hence, if G^ = gH£ u , 

d*G = *j mag . (3.4) 
Dirac defines the electromagnetic field by 

F = da + *G (3.5) 
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and gets the desired modified Bianchi identity dF = *j mag - Note that we have used the 
lowercase for the electromagnetic potential here. This is to distinguish it from the 
potentials Af in our formalism ( see Eq. (12.11) ). In particular, A\ in our formulation is not 
equal to tij, which arises in other two-potential formulations to be discussed below. 

In the Dirac formulation, the theory has an extra gauge symmetry associated with 
the freedom of arbitrarily choosing the position of the strings while keeping its boundary 
(worldline of monopole) fixed. To see this, consider the displacement of a string defined 
by 

x" = w ii (t,(t,\). (3.6) 

where the initial string worldsheet £ is at A = and the final, £', at A = 1. The boundary 
of the 3-dimensional surface T defined by (13.61) is AS = £ — £'. Hence, if 

AHg = Hg - (3.7) 

then from fO) . 

*AG = d*K, (3.8) 

where K a ^ = gHf 1 '. Therefore, we see that the electromagnetic field F in (13.51) is 
invariant under the displacement of the string if, while moving the string, we also vary a 
by 

Aa = -*K. (3.9) 

The Dirac action, which depends on the string only through F^ u is invariant under this 
gauge symmetry, up to the anomaly that gives rise to the quantization condition, which 
will be explained in more detail below when discussing the double potential formalism. 

In a manifestly duality invariant theory, magnetic and electric charges are treated on 
the same footing. In general one considers n dynamical dyons with magnetic and electric 
charges q% = (g n , e n ). The current is then defined as 



where the sum in n is over the worldlines T n of every dyon of charge [parameterized 
by z£(t) with an arbitrary parameter r]. For the Dirac strings attached to them, we define 

GW{z) =J2l a nHZ(x) = I ~ Vn)dyZ A dy v n , (3.11) 

n n 

where £ n is the worldsheet of the Dirac string whose boundary is T n [parameterized by 
y%(r, cr) with arbitrary parameters r and a]. The analogs of Eqs. (13.41) and (|3.8I) in this 
case are 

d*G a = Y, *AG a = d*K a , (3.12) 
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where 

K aa^ = ^ q a H *(h ? (3 _ 13) 

n 

and T n is the surface defined by the displacement of the string attached to the dyon q%. 

When splitting space and time, as in the different manifestly duality invariant formu- 
lations, it is convenient to also split the space and time components of the string currents, 
defining 



, o, 



2 f 



/~ia jk */~ia nai /~ia Oi ijk */~ia /o i a\ 

a i = ~ e ijkLr = G 0i , p — G = -e J G jk (3.14) 



Explicitly, 



a" 



^2<ln S 4 (x- y n )\dy n x Ady n , (3.15) 
P a = X>n / S\x-y n )dy° n Ady n , (3.16) 

where (dy n x Ady n )i = e^dy 1 A dyi and the first identity in (13.121) becomes 

f ■ f3 a = j a0 , V x a a -d (3 a = f. (3.17) 

It is also convenient to work with the dual of K a a/?7 , the one-form v a a . In terms of it, we 
may derive the way the vectors a a and (3 a in (13.141) transform under displacement of the 
strings. Using the second identity in (13.121) . 

Aft? = *AG a 0l = (dv a ) 0l = d v?-divl (3.18) 

Ap ai = \e ijk *AG% = \e i]k (dv a ) jk = (V x (3.19) 
2 J 2 

For dynamical dyons, the action must be supplemented with the kinetic term 

hK] = - [ yj-dz£dz nii . (3.20) 

The total action I' that includes the dynamics of the dyons and produces the correct 
Lorentz force law is 



I'[AlC a ,y^=I M + Ii + h+ 



+ - j d '.r £„,, 



2VC a a b - (3 a a b - /3 a V~ 2 V x d (3 b 



(3.21) 



The constraints (12.71) and the electromagnetic Eq. (12.91) are clearly unchanged, for the 
extra piece in the action does not depend on A" The equations obtained from the variation 
of C a are modified with respect to the result of (12.81) to 

V 2 C7 a = e ab (V ■ [dvJP + ft 6 ] - V 2 ^) . (3.22) 
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For later use we note that, applying Vx to (12.91 ) and using (|3.22l) . together with the 
boundary condition that B a falls off at least as fast as r _1 at infinity, 

B a « e ab (d A b -VA b + a b + V- 2 V xdo^y (3.23) 

As a side remark, we also note that the definitions 

K» = d,At-d u A^ + *g; u , *g*=e l]k d k C a , (3.24) 

*Q«. = ft «+V- 2 (Vx9 ft. (3.25) 

are such that 

B a i = -eijkF{ k = *F a0i . (3.26) 
Furthermore, they allow us to write the on- shell Eqs. (13.231) as 

B a i ~ eafeFgj = —-^eij k e ab *F jkb , (3.27) 
while the equations of motion for A a ^ take the covariant form 

d v F>? = e ab j b » t 0„*Ff = -#. (3.28) 

In the case of a single dyon, one assumes without loss of generality that the string 
terms in the last line of (13.211) are absent. Indeed, in this case one can perform a du- 
ality rotation so that, say, the magnetic charge vanishes. The string terms then reduce 
to — J d 4 xVC l a 2 and can be dropped because they only affect the irrelevant auxiliary 
equation used to determine Aq. This justifies a posteriori the coupling to the sources 
considered in the first section. 



3.2 Lorentz force law and veto 

We still need to vary y^, in action (13.211 ) in order to derive the Lorentz force law. We 
will see below that in order to obtain it, we need to impose the so called "Dirac veto." 
This demand was introduced by Dirac in his original treatment of magnetic monopoles 
lfT2l to obtain the desired classical equations. It consists of the requirement that no electric 
charge can touch a Dirac string. At the quantum level, Dirac showed that the veto modifies 
the topology of phase space, giving rise to his celebrated quantization condition. In our 
formalism the Dirac veto is required as well, as we show below. The difference resides 
in that here we will need to ask that no dyon can touch the string of any other dyon. This 
generalized version of the Dirac veto was also used in J9]|. 

Variations of Ij with respect to z£ give 

SJi = e ^n I {(9 Q A a - VA a ) ■ (5z° n dz n - 5z n dz° n )+ 

+ (V x A a ) ■ (5z n x dzS) ■ ( 3 - 29 ) 
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Before varying Im with respect to y£, we establish the following identities. For all 
smooth vector fields V a , W a one has 



J d 4 xV b 5 y a b = J2<lr, 



1 



V b ■ (5z n x dz n ) + I V ■ V b 5y n ■ - (dy n x Ady n )- 



- d V b ■ [5y n x {dyl A dy n ) - 5yl\ 



x Ady n ] 



(3.30) 



and 



J d'xb y ^W a = Y,<fn 



(8z° n dz n - 5z n dz° n ) ■ W a + 



+ / (Sy n x (dy° n A dj/ n ) - 8y Q n \ (dy n x Ady n ) • (V x W a )j 



. (3.31) 



The variation of I M with respect to y^ may be computed by specializing for the fields 

1 



V b 

Wa 



■e ab (2VC a - (3 a ) 



l 



--e ab (a b + 2V- 2 Vxd (3 



(3.32) 
(3.33) 



Combining all terms, 

SJl + Syhl = ^ l ((^a- e ab (d Q A b ~ V^)) ■ (6z° n dz n - 6z n dz° n ) + 

+ (Va ~ e afe (V x A 1 )) • (5z n x dz, 

+ (V x W a - d V a ) ■ Uffn x (dy Q n A dy n ) - 6y°~ {dy n x Ady„) 



V • V a 8y n ■ - (dy n x Ady n ) + 



(3.34) 



Now, taking the divergence of (13.321) and making use of the first identity in (13.171) and the 
constraints (12.71) one gets, 

1 



V • V a = - e ab j 



60 



(3.35) 



It follows that the second term in (|3.34l) vanishes provided the string attached to dyon n 
does not cross any other dyon (Dirac veto). This is due to the fact that the Dirac veto 
ensures that j 7 " 1 = on the worldsheet of the strings. Similarly, from (13.321) . (13.331) and 
the identities (13.171) it is straightforward to show that 



V xW a -d Q V a = -\e ab f. 



(3.36) 
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Hence, the last term in (|3.34l) also vanishes on account of Dirac's veto. The string piece 
in the first and second terms of (13.341) again vanish because of the veto. This may be seen 
from the fact that (13.321) may be written as, 

V a - e ab (V xA b ) = -e ab B b + \e ab f3\ (3.37) 

and therefore, due to the veto, the integral on the worldline of a dyon only sees the first 
term. In the same way, using (13.231) . 

W a - eatidoA" - VA b ) = -B a + \e ah a\ (3.38) 

and the second term vanishes on the worldline of a dyon. 

Combining the remaining terms with those from the variation of I k , extremization of 
the total action now implies the Lorentz force law 

m n j- f j2L = ) = q a n B a (z n )-f^, (3.39) 

CLT \ / dz nu dzX / GST 





dz£ 


V dr 


dr 


dz„ 




dr 




/ dznu, 


dz% 


V dr 


dr 



m n - , dT = q a n B a (z n )--^ + -^xB a (z n )e ab ql (3.40) 

dr \ / dz nu dzii / dr ds 



as it should. 



3.3 Equivalence with Dirac's co variant formulation and quantiza- 
tion condition 

We end this section by showing that the theory presented above is equivalent to Dirac's 
theory. This shows that the theory with dyons is Lorentz invariant. We will actually show 
that our action (13.211) is equivalent to an action found in which, in turn, has been shown 
in ifTOll to be equivalent to a generalization of Dirac's covariant formulation allowing for 
dyons. 

Explicitly, this action reads 



I[a a , y§ = \ J d 4 x [ea^idoa* + a b ) -b a -b a + t ab d a ■ f 



+ h, (3.41) 



where 

b a = v x a a + /3 a . (3.42) 

This formulation makes use of Dirac strings in the same way as our formulation does. 
That is, each dyon q n is attached to a string parameterized by y£(r, a). The quantities a a 
and (3 a appearing in (13.411 ) are the same ones as defined in Eqs. (13.151 ), (13.161) above and 
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satisfy the identities (13.171) . Note that on account of these identities, the longitudinal part 
of a drops out of this action principle. The field b a is the magnetic/electric field appearing 
in Maxwell's equations. It must, therefore, be the same as our B a . 

— * 

In this formulation, Gauss' law appears as an identity on taking the divergence of b a 
in (13.421) . The field V x a a is transversal but has a stringlike singularity which is removed 
by (3 a . In our formulation B a has two nonsingular pieces, namely, the transverse and 
longitudinal components of it. To show equivalence, we decompose (3 a accordingly so 
that 

fe a = Vxa a + /3 a = Vx(a a - V~ 2 V x f3 a ) + VV~ 2 V • (3 a . (3.43) 

From the constraints (|2.7I) and the first identity in (|3.17l) . the longitudinal piece is precisely 
VC7 a . We are then lead to the following identifications: 

A a = a a ~ V~ 2 V x /3 a , (3.44) 
C a = V~ 2 V-/3 a . (3.45) 

(The first equation is true up to an irrelevant longitudinal field, in order for B a = b a to 
hold.) 

To establish equivalence is now straightforward. We start with the action (13.211) of 
our formulation. Assuming vanishing boundary conditions at spatial infinity, (Aq, C a ) 
are auxiliary fields because their equations of motions (12.71) and (13.221) can be used to 
algebraically determine them in terms of the other fields. We can thus solve for them in 
action (|3.21l) . Then we use (13.441) to write A a in terms of a a , and after a bit of algebra 
involving the identities (13.171) we get precisely action (13.411) . 

In the double potential formulation of flU, i- e -> for action (|3.41l) . the symmetry cor- 
responding to shifts in the string is realized by transforming a a and (3 a in (13.181) . (13.191) . 
with v o = and 

Aa a = (3.46) 

Let us do that for the case in which there are only two dyons, q a , q a . We will only change 
the position of the string attached to q a . Varying action (13.411) and using the identity 

V • b a = j a0 we get 

5I=\tab j d*xj*»vl (3.47) 

This is zero unless the worldline of dyon q a crosses the 3-dimensional manifold T swept 
by the string attached to q a . In that case the variation is 

5I= l -e ab q a q b . (3.48) 

This will not affect the quantum mechanical system if the variation is proportional to 
2nhn, for some integer n. This leads us to the Dirac-Schwinger-Zwanziger quantization 
condition 



eg — eg = 2irnh, 



(3.49) 
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up to a factor of 1/2. This factor is removed by a careful analysis of the topology of the 
system. We will not discuss this here. More details can be found in [|9l ITOTl. 

Finally, we study what happens in our formulation. First, let us compute how the field 
A a transforms under the movement of the string. From Eq. (13.441) . (13.191 ) and (13.461) we 
get (EH, 

AA a = -v a - V~ 2 V xVx)f = -V (V~ 2 V • u°) . (3.50) 

If we now take 

AA a = -w a = -do (V 2 V • + d (V 2 V • - vg, (3.51) 

the variation defined by (|3.50l) and the first term of (13.511) is a gauge transformation of 
the form (12.61) which leaves action I' in (|3.21l) invariant. We thus only need to compute 
the variation under the movement of the strings and the second part of (|3.51l) . Using 
identities (13.171) one obtains precisely the same result as in the previous case, namely, 
the right-hand side of Eq. (13.471) . The argument leading to the quantization condition is 
therefore the same. 

4 Extended double potential formalism in curved space 

We generalize the first order action to curved spacetimes and discuss the canonical and 
gauge structure of the theory, including diffeomorphism invariance. In particular, we 
show that the standard algebra of surface deformations of the purely gravitational case 
now involves both Gauss-type constraints with structure functions depending on electric 
and magnetic fields. We proceed to the equations of motion deriving from the general- 
ized action principle and show that they are equivalent to the covariant Einstein-Maxwell 
equations. We then show how the string singularity of the Reissner-Nordstr0m dyonic 
black hole solution gets resolved in our formalism. We compute the electric and magnetic 
surface integrals following the Regge-Teitelboim approach, discuss how they appear in 
a geometric derivation of the first law and in the Euclidean approach to black hole ther- 
modynamics. Finally, we apply these results to the resolved Reissner-Nordstr0m black 
hole. 

4.1 Action and canonical structure 

The first order action Im can be generalized to curved spacetimes. We consider a globally 
hyperbolic spacetime, foliated by a spacelike family of hypersurfaces, each labeled by 
the value of a timelike coordinate t. The induced metric on each surface is gij(t). We 
follow the conventions of MTW [fT5l . chapter 21, where spatial indices are lowered and 
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raised with the 3-metric and g is its determinant. We denote by e^k the completely 
antisymmetric symbol, which differs from the [ijk] notation used in MTW. 

Adapting the results derived in (181171 l8ll6ll and defining B ai = e^ k djA a k + y/g&C a , 
we get the following manifestly duality invariant action in the absence of sources: 

I M [A^ C\ 9ij , N, N l ] = ^J d 4 x [(B ai + ^C a )e ab (d Al - d^)- 

-^-BiBl-e^N^B^}, (4.1) 

where iV = (— Wg 00 ) 1 / 2 and N l = ^'goi are the lapse and shift functions and ^g^ v 
is the 4-dimensional metric. 

We are interested in solutions to the equations of motion derived from I = Iadm+Im, 
where Iadm is the first order action for pure general relativity. Introducing the collective 
notation z A = (g^, Af, C a ) for the different fields in our system, this action principle 
takes the form 



I[z,u} = J d A x [cLA(z)d z A — w a 7a], 



(4.2) 



a A {z)d z A = ^-d o9ij - ^d A t + y^^d Z t . (4.3) 

1D7T 47T 47T 

The constraints 7 Q = (H±,Hi,G a ) are associated with the Lagrange multipliers u a = 
(N, N\ A$) and given by9 

H± = -^-{Hl DM + HT), Hi = -^-(Hf DM + HT at ), Qa = ^e ab diB b \ (4.4) 
lo7r " lo7r Air 

where Hf DM , Hf DM are given in(I51|Bl and 

n rnat = ^gig<y ? n ™t = 2 e ab e ijk B a3 B bk . (4.5) 

The first two sets of constraints in (14.41) above are the gravitational Hamiltonian and 
momentum constraints, while the last set are the two electromagnetic Gauss constraints. 

In order to disentangle the canonical structure we begin by writing this action as 



l -Jd A x[- £%A + y/g&CdoZi - Aod.S 1 + ZodiB'- 

(4.6) 



I M = I d x 



1 ^(8 l 8 l . + B l B i ) + e l3k N l 8' J B k 



'Note the misprint in eq. (21.116) of 02), where there should be no lapse function on the right-hand 
side. 
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where E % = e ljk djZ k + Jg&'Y and B 1 = e %jk djA k + y/gd l C. We assume here and below 
that every 3-vector admits a unique orthogonal, spatially covariant decomposition (see 
e.g. [17]) X i = X Ti + X Li , where 

X Li = &M, X Ti = —e ijk diL k , (4.7) 

for some M, L k . In terms of the inverse of the spatially covariant Laplacian V~ 2 and the 
spatially covariant derivative V*, we have 

M = S7~ 2 VjX j , X Ti = X i - d { M. (4.8) 

A vector is transverse if its divergence vanishes and longitudinal if its curl vanishes, 

dii^/gX 1 ) = => X i = X T \ e ijk djX k = => X i = X Li . (4.9) 
We then have 

J SxyfgX'g^ = J d 3 x^(X Li gij Y L i + X Ti gij Y T i). 

Using such a decomposition for A\, A* = diM a + Af , M a = (M A , M z ), the kinetic 
term becomes 

fd i xa A (z)d ri z A = J d l x([^ + ^ (Z% V9d>C - Aljgd,Y)] %g tj - 



47T An An 

where D ijH = ^=(g ik gji + g u g jk - g^gu) is the DeWitt supermetric. Note that D l]kl is 
not the inverse DeWitt supermetric, but the result of raising all the indices with the metric. 
Let us define 

vr a = (di(- s /gd i C), di( y fgd i Y)) = (-n z ,n A ), (4.11) 
as new independent variables, so that C a = V~ 2 ^=. We also define 



n 13 = n l] 



+ A^D^ kl e ab Af ^dtC\ (4.12) 



The independent phase space variables are thus (g^-, Af , M A , M z , n l \ Zj, n A , t^z) m 
terms of which the canonically conjugate pairs are 

~kl c ijkf) 7T 

^•k>' ^v^po, ( M -S»' <4 - 13) 

In particular, 

{B ai (x),B b \y)} = Ane^ k e ab d x k 5 3 (x,y), {M a (x) , n b (y)} = Ane ab 5 3 (x,y). (AAA) 
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4.2 Gauge structure 

Before turning to the equations of motions and their solutions, let us discuss the gauge 
structure of the theory. We want to show that the constraints j a are first class. Defining 
e a = -^ a ) with e a vanishing at the boundary and r[e] = / d?x^j a e a , this means 

that 

{r[e 1 ],r[e 2 ]} = r[[e 1 ,e 2 ]] > (4.15) 

for a suitably defined [ei, e 2 ]. In this case, the gauge transformations leaving action (14.21 ) 
invariant are given by 

5 e z A = {z A ,T[e}}, 5u a = d e a + [e,u] a . (4.16) 

In order to compute these brackets it is useful to go to the Darboux coordinates iden- 
tified in the previous subsection in terms of which the Gauss constraints become 

G a = ^e ab n b . (4.17) 

47T 

Now one should do the change of coordinates in TC±,TCu i.e., perform the replacement 
n ij _ _ /±^~gD l i kl t ah A°^ \j~gdiC h '. Since the additional terms that are generated in this 
way are proportional to C a and thus vanish on the constraint surface defined by Q a they 
can safely be discarded in the source-free situation. In the following, we will drop the 
tilde on 7r u . 

In particular, we have 

5 e ir a = 0, (4.18) 

5 e A° = d i \ a T -^e ab Btf ± -e ijk ¥B ak , (4.19) 
v 9 

Se9ij = V i £ J + V£ i + 2D tjkl 7r kl £ ± , (4.20) 

where 

X a T = X a - e^V- 2 ^^ 1 ) + -iv-^^^). (4.21) 

Equation (14.181) implies that n a are constants of motion, which is consistent with the 
longitudinal part of the source-free Maxwell equations. From (14.191 ) we also find 

5 e B ia = -e ijk dj{^-e ab B bk ^) - d^C) + dj{B ai ^). (4.22) 
yd 

Infinitesimal diffeomorphisms along i]^ are recovered by using = Nrf, & = g^rf 1 . 
Indeed, with this choice of parameters, C^g^ ~ S^g^ v . This can be seen for instance on 
(14.201) by using the (auxiliary) equations of motion for rr^ together with the definitions of 
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lapse N and shift N l in terms of the 4-metric g^ u . In other words, diffeomorphism invari- 
ance in the Hamiltonian framework is implemented through the gauge transformations 
generated by H[£\ = f d 3 x (H±£, ± + TiiC)- I n particular, using as gauge parameters the 
Lagrange multipliers, e a = u a , amounts to performing an infinitesimal time-translation 
on account of the Hamiltonian equations of motion. 

Let us end this discussion by determining [ei, e 2 ]. The constraints Q a , have vanishing 
Poisson brackets among themselves and with all other constraints because the B ai do not 
depend on M A , M z . It follows that [A, e 2 ] a = and also, from (14.191 ) and (14.161) that 
the associated gauge transformations 5\ generated by G[X] = f d 3 x Q a ^ a are the double 
electromagnetic gauge transformations of (12.61) . while all other variables are left invariant. 

We still have to compute {#[£], H[i]}}. We notice first of all that the purely gravita- 
tional part satisfies the algebra of surface deformations [fT8l[T9l . 

{H ADM [$, H ADM [ V }} = H ADM [% V } SD ], (4.23) 
[£,ri]sD = CdW L -V i di£ x , (4.24) 
v]sd = tfK^drf - V ± d^ ± ) + t'drf - rf dtf. (4.25) 

From (H32l) and (l4~20l) . we find that 

{H ADM [£}, H mat [7]]} -(Z++ v ) + {H mat [^},H mat [r]]} = 

= H mat [[^ V ] SD ] + G[[^r ] ] B ], (4.26) 

where 

f ac B ■ 

vTb = B a %jk{ j V k - -^(ZW - V X C)- (4-27) 

v 9 

Combining with (|4.23l) . we finally get 

{H[Z],H[tj]} = H[[t,7j\ SD ] + G[[£,rj\ B ]. (4.28) 

According to [|20l . such a constraint algebra provides the integrability conditions that 
guarantee that "the evolution of a three geometry can be viewed as the deformation of a 
three-dimensional cut in a four-dimensional space-time". 



4.3 Derivation of the Poisson algebra of Poincare generators in flat 
spacetime 

In this subsection, we derive the Poisson algebra of the Poincare generators in flat space- 
time as given in (12.231) by restricting the results of the previous subsection to flat space- 
time. 
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We thus assume in this subsection that N = 1, N l = 0,gij = Greek indices 

take values from to 3 with jx = (±, i). Indices are lowered and raised with r/^ u = 

diag(— 1, 1, 1, 1) and its inverse. Let ui^ = —Lo vtl . In this case, the Lie algebra of vector 

_ ^ _ _ g 

fields ^(lo, a) = (— u^x 1 + a ^)g~^i with bracket the surface-deformation bracket (14.241) 

and (14.251) forms a representation of the Poincare algebra [EH, 

[£(£i,ai),£(£3 2 ,a 2 )]sD = £(Pi> ^2], ^i«2 -^201). (4.29) 

It then follows from (14.281) that this is also the case for the canonical generators ?i[£(u, a)} 
equipped with the Poisson bracket, when one considers the restriction to the constraint 
surface defined by Q a = 0. 

Comparing with Sec. 12.41 we find that 

~uvJ^-a^ = W[£(£,a)] 

Ufu, = iTTUfja,, a = 4na , a { = A-n(di + UoiX ). (4.30) 

This concludes the proof that the generators defined in (12.221) form a representation of the 
Poincare algebra and the algebra in (12.241) is a direct consequence of (14.281) . 



4.4 Equations of motion with sources and comparison to covariant 
formalism 



The standard Einstein-Maxwell equations, now in the presence of external, magnetic and 
electric conserved current densities j af *, c^j aAt = given by (13.101) with associated string 
terms defined in (|3.1 II) . derive from extremizing the action I geom 



I' M , where 



L geom [tifii/ 



(4.31) 



and 



1 

47T 



(fx 



- 2 ( d l* a " 



1 



1 

+ - 



(4.32) 



4^ 

To make connection with the ADM formalism, we have followed [fT6ll and introduced the 
auxiliary tensor densities T^ h '. On the one hand, one can solve the equations of motion for 
^algebraically,^ = a/ -( 4 > gg tMa g u(3 (d^a u -d I/ a fM + *G fMU ). When substituted into the 
action, one recovers the Einstein-Maxwell theory with Dirac strings. On the other hand, 
one can introduce F 0t = £\ DM , eliminate the auxiliary T % i and use the decomposition of 
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the 4-metric into the 3-metric, lapse and shift to find the standard Hamiltonian form 

AA^ADM-, <V 9ij, N , N \ yl = ^ / d 4 X - £ ADM {d ai + a i) ~ a O d i^ADM- 



^ ( P { cADM _| joi joADM \ , Art pi »/s , 

l^ADA/^i ' °ADM°i I ~r e ijfc JV O adm O adm + <2 M J 



(4.33) 



2v^ 

where B l ADM = e^ k dja k + (3 k . The equations of motion for £ ADM read 

^adm = ^-9 ij (dja - dodj - a, - e jkl N k B l ADM ). (4.34) 

They determine £ ADM in terms of the other variables and the sources. Similarily, in the 
gravitational sector, the equations of motion following from varying 7r u are auxiliary in 
the sense that they can be solved algebraically for 7r iJ in terms of the other variables. After 
this has been done, the constraints and the equations of motion following from variation of 
gij are equivalent to the covariant Einstein-Maxwell equations with Dirac strings. Hence, 
every solution g^,a^ to the covariant equations of motions is a solution to the ADM 
equations of motion with electric and magnetic fields £ ADM , B l ADM and momenta it 1 ! 
determined in terms of g^ u , a M . Conversely, every solution of the ADM equations of 
motion gives a solution to the covariant equations of motion. 

Alternatively, one can multiply (14.341) by N. The longitudinal part of this equation is 
solved uniquely for ao, while the transverse part gives, after using e %ik dja k — do(3 % = k\ 

doB ADM = -e^(^ DM ) - d 3 (N*B ADM ) + d {WB\ DM ) - k\ (4.35) 

Finally, Maxwell's equations for Oj are 

do£ AD M = e^%(^,Bt DM ) - <9,(iV^ DM ) + d^£ ADM ) - j\ (4.36) 

As a side remark, note that the constraint algebra in the absence of sources and strings 
in the standard ADM approach to Einstein-Maxwell theory can be directly rederived from 
our result (14.281) and is given by 



{H[£],H[7j\} = H[[£, V } SD ] + G ADM [[^ V ] ADM ], (4.37) 



where 



Gadm[X] = J d 3 xdi£ ADM \, 



cADM 

%v\adm = SUaV-V^"^' (438) 

Indeed, the algebra rests only on the constraints and the transformation properties (|4.20l) . 
(14.141) . Provided that £ ADM = £ % and B\ DM = B\ which we always assume in the fol- 
lowing, these are the same in both descriptions, except that the constraint Q 2 = —diB ADM 
is absent because diB ADM vanishes identically in the absence of sources and strings. 
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In the presence of sources, (13.411) generalizes readily to curved space, where a mani- 
festly duality invariant action principle is defined by Iadm + Im, with 

I M [al 9ij , N, N\ y^} = -L / d 4 x \(b m e ab (d a1 + a\) - (4.39) 

- e ab e l3k N^b hk + e ab a<Zf] , (4.40) 

and 6? = e l:jk djat + (3 a% . Indeed, equivalence of the associated equations of motion to the 
ADM/co variant ones is obvious when b 1% = B\ DM , b 2t = S l ADM since the gravitational 
equations of motion are unaffected while those for a" read 

d b ai = -e^dji—e^hk) - dj(N l b aj ) + dj(N 3 b ai ) - f\ (4.41) 

and coincide with the relevant Eqs. (|4.35l) - (14.36l) . 

With the longitudinal electric and magnetic fields produced by the potentials C a , the 
appropriate action principle is Iadm + Im + Ij, where Im is defined in (14.11) and 

Ij[Al,C a ,y»] = ±-J d*xe ab (A^ + ^d*C a at- 

-l^al + l^d^l). (4.42) 

Here 7" is the potential for the transverse part of (3 a% , (3 aT% = e^ k dj^. 
In this case the equations of motion for A" are given by 

d B ai = -e^dji—e^Bbk) - dj(N l B aj ) + dj{N j B ai ) - j ai . (4.43) 

They are the correct matter field equations provided that B m = b ai are the magnetic and 
electric fields. This implies on the one hand e^djA^ = e^ k dja1 + (3 aT% , and in turn 
At = a k + 7fc> U P t0 an irrelevant longitudinal part, and ^Jgd l C a = (3 aLt on the other 
hand. Again, the equations of motion for Aq, C a , 

di(Jg#C a ) =j a0 , (4.44) 



di 



0, (4.45) 



are auxiliary because they can be used to solve these fields in terms of the others. This 
can be done in the action principle and gives back (14.391) . 

In conclusion, if the lapse N is nonvanishing and the covariant decomposition of 
spatial vectors into longitudinal and transverse components is unique then there is a one- 
to-one and onto correspondence between solutions of the covariant/ ADM equations of 
motion and solutions to the equations of motion deriving from Iadm + Im + I J- 
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In the case of a single dyon, one can again drop all string terms in Ij, which then 
simplifies to 

Ij^ri = ^ J d 4 xe ab A*j b ». (4.46) 

All relevant matter equations of motion are correct in this case, but one has to face the 
fact that the metric dependence in the longitudinal part of B m implies an additional term 
in the equations of motion associated with variations of g k i, 

8(1 ADM + Im) _ 5(1 ADM + hi) _ \f~9~ jjjjkl q (ja-s^ 

Sgki 8g ki j ai, 

K = ^-B\ - e ab g il (d A b - d t A b + e lkm N k B bm ) . (4.48) 

Again, one can use a duality rotation to make the magnetic charge vanish in which case 
the equations of motion imply C l = 0. We thus only need to consider X\. But in the 
purely electric case a« = and = so that X % 2 vanishes on account of the matter 
equation of motion (14.341) . 



4.5 String-singularity free dyonic black holes 

Consider now the case of a dyonic Reissner-Nordstr0m solution with charge Q a . The 
dyon defined by (ll.ll) - (ll.2l) is a solution outside of the location of the dyon (at r = 
in our coordinate system) with a Dirac-string singularity to the equations deriving from 
Igeom + I'm gi ven in (14.311) and (14.321) . It is thus also a solution to the equations of motion 
derived from Iadm + I'm f° r which 

4dm = KQ sin 9, B\ DM = 5).P sin 9. (4.49) 

In the simplified duality invariant formulation defined by I ad m + hi + Ij, with Ij 
given in (14.461) . we have to determine the vector and scalar potentials giving rise to B ai = 
5 l r .Q a sin 9, where Q 1 = P,Q 2 = Q with the metric given by (11.11) . This is easily seen to 
be the case for A a% = and 

dr' 

C a = -Q a ■ 



r' 2 N(r') 



Q a ln r(M-yWQ~f) 



^fQK[ b Mr - Q-Q c - ^Q d Q d (r 2 - 2Mr + Q e Q e ) 
-— + 0(r~ 2 ). (4.50) 



r 

In the gauge where the scalar potentials vanish at infinity, it is then straightforward to see 
that all matter equations of motions are solved by 

At = (4.51) 
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and one can directly check that in this case X l a = 0. In conclusion, in the new formulation, 
the Reissner-Nordstr0m dyon is described by the metric (11.11) and the potentials (14.501) . 
(14.511) . The string singularity of the standard approach has thus been resolved in the new 
formulation. 

In the gauge where the scalar potentials vanish at infinity, let us define 

<j) = -A), V = Zo- (4-52) 

with (j>Hi ipH denoting these quantities evaluated at the horizon, in agreement with (11.61) . 
For the resolved Reissner-Nordstr0m dyon, this gives 

Q P 

In the Euclidean methods discussed below, it is useful to choose a gauge where the scalar 
potentials vanish on the horizon. In this case, 

A a = -e ab Q b (- - — ), (4.54) 
r r + 

A = -^ + <fe, Z = j-iP H . (4.55) 

4.6 Surface charges 

The Regge-Teitelboim analysis [f2TTl allows one to derive the correct variational princi- 
ple in the presence of nonvanishing surface charges at infinity. Consider an arbitrary 
gauge transformation generator r[e] = J d z x^ a e a with e a not necessarily vanishing at 
the boundary (a condition we demanded in Section 4.2 above). The variation of this 
generator under a change of phase space variables may be written 

5 z T[e] = J d 3 x5 z ( la e a ) = J d 3 x (<^ A %P - , (4.56) 

where 5/5z A is the Euler-Lagrange derivative. The second piece is a boundary term and 
arises from integration by parts. The expression k\ [z A , Sz A ] depends on the phase space 
variables and linearly on their variations and the gauge parameters. 

For the simplest application, consider phase space variables z A that satisfy the con- 
straints and variations Sz A obeying the linearized constraints. In this case, the left hand 
side of (I4.56I) vanishes. Suppose then that the associated solution z%, it" to the evolution 
equations is time independent, d z A = 0. In particular, this means that the associated 
vector field = 5$ is the timelike Killing vector field of the metric g^ u . In this case, the 
evolution equations following from (I4.2I) . 

gf^-a^ = fe^>, (4.57) 
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imply that the first term on the right hand side of (14.561) vanishes as well. We thus find 

Using Stokes' theorem, it follows that the integral over a sphere at radius r and fixed time 
t does not depend on the radius r, 

I d 3 Xi ki s [zf,5zf]= I d 3 x t kl[zf,5zf]. (4.59) 

Here d 3 Xi = \e^ k dx^ A dx k . The explicit expression for k\[z A ; 5z A ] can be easily worked 
out by integrations by parts. It is defined up to the divergence of an arbitrary superpoten- 
tial, djt™, which does not play any role for our purpose. It splits into a standard purely 
gravitational part and a matter part, 

k\[z A - 5z A ] = kr v '%j, * lJ ; S9ij, S7r ij ] + fc^Vl 5z A \. (4.60) 
The former has been derived in [|2T1l and reads 



J^grav,i 1 

e 16tt 



& k \ev k 5 9lj - d k ^5 gij ) + 2£ k 5n ki + (2£*7r* - £V fe )<% fc 
G ljki = \^g{g lk g jl + g a g lk - 2g lj g ki ), 



(4.61) 



where G ijkl is the inverse of the DeWitt supermetric, D ijkl G klmn = § + 5f5f). The 
matter part now involves, besides the electric contributions, the sought for magnetic ones: 

k vaat,i = (^_ e ijk B a 5Aak _ ^gai^ + e^B™ - CB ak )8A\- 

- e aby /gg% k eB ak 6C b + e ab { s fgd i X a 5C h - \ a 6B bU )) . (4.62) 



In the sequel, we are interested in asymptotically flat gravitational field configurations car- 
rying finite charges associated with energy momentum. We will not need to consider the 
more general boundary conditions guaranteeing finite charges associated with rotations 
or boosts. The appropriate fall-off conditions on the gravitational variables and lapse and 
shift have been discussed in detail in 112111 . 

g rr = 1 + 0{r~ l ), g ee = r 2 + 0(r), g u = r 2 sin 2 6 + O(r), (4.63) 

9re = O(r ) = g r +, ge<f> = 0(r), (4.64) 
n rr = O ( r 0) ; n ee = ^ = n H = ^ = Q ^ = ^ (4 65) 

N = 1 + 0(r~ l ), = 0{r- 2 ) = N 9 , N r = 0(r -1 ). (4.66) 
For the matter variables, we assume 

A* = 0(r~ l ), Ag = 0(r°) = A^, C a = 0(r _1 ), A a = k a + 0(r~ l ). (4.67) 
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In particular, these fall-off conditions include the background solutions z, u described by 

g rr = 1, gee = r 2 , g H = r 2 sin 2 9, (4.68) 
N = 1, N 4, = = N e = N r , Al = k a , (4.69) 

and all other variables vanishing. For later use we introduce the additional notation 

k 1 = C , k 2 = -^ c . (4.70) 

In order to allow configurations satisfying the fall-off conditions to be extrema of the 
variational principle, action (|4.2I) needs to be supplemented by the addition of a suitable 
surface term at the boundary at infinity, i.e., the surface r, t constant with r — > oo, 

I T [z,u] = J d A x[a A (z)d z A - u a la \ - Q u [z]. (4.71) 

The surface term Q u [z] is determined by the requirement that, under variations of the 
fields z A satisfying the fall-off conditions, its variation 5 Z Q U should precisely cancel the 
spatial boundary term arising when deriving the Hamiltonian equations of motion, i.e., 
the term due to the right hand side of (14.561 ), 

6Mz] = <f iPxihifoSz]. (4.72) 



For the purely gravitational part, this problem was solved in [21 J, the appropriate bound- 
ary term being the ADM mass: 

d 3 Xi fc£ rw ' < [z, 5z] = <f d 3 Xi k 9 u rav,t [z, 5z] =8 z <f d 3 Xi k 9 Q rav,i % z - z], (4.73) 

S°° J S°° J S 00 

so that 

Q*T [g,7r]=<f d 3 Xl k 9 u rav>i [z } z-z]=M, (4.74) 
M = I d 3 x t V§(9 lk g ji ~ 9 lj 9 ki )D k (9ij ~ 9ij), (4-75) 

where the covariant derivative is taken with respect to the flat background metric g~y. 
For the matter part, the boundary conditions imply also that 

d 3 Xi C aM [z, Sz] = I d 3 Xl k™ at,i [z, Sz] =5 z (f d 3 Xi k™ at,i [z, z-z\. (4.76) 
Js°° Js°° 

In particular, the boundary conditions (I4.63I) - (I4.66I) are such that, when (^" L ,^*, A a ) are 
replaced by (N, N\ A%), the contributions proportional to N,N l from the matter part 
(14.621) vanish. There is thus no correction to the ADM mass for the adopted boundary 
conditions. This will not remain true for more general boundary conditions where the 
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matter part (14.621) can contribute both to the ADM energy momentum and the Lorentz 
generators. For the boundary conditions at hand, only the last term survives and combines 
into magnetic and electric charge Q a = (V, Q), as expected, 



Q™ at [g, C] = j> d 3 x t z-z} = -k a e ab Q\ (4.77) 



Qb = J- ( f ) d 3 Xl B biL . (4.78) 

47T 



In other words, off-shell, the correct Hamiltonian for the boundary condition under 
considerations is 



H = J d 6 x (H±N + KlN 1 ) + M, (4.79) 
while electric and magnetic charges are given by 

Q = -j c j d 3 x (SiAo) + Q, P = ~J d 3 x (G 2 Z ) + V. (4.80) 
These obervables commute in the Poisson bracket, 

{H, Q} = = {H, P} = {Q, P}, (4.81) 
and the total action (14.711) can be written as 

I T [z, u] = f dt(J d 3 x a A (z)d Q z A - (H - <p c Q - ^P)^ . (4.82) 
4.7 First law 

For the resolved Reissner-Nordstr0m dyon z, u given by (11.11) . (14.501) . (14.551) . the first law 
of thermodynamics can now be derived as a consequence of using identity (14.591) between 
infinity r\ — > oo and the outer horizon r 2 — r + , 

I d 3 x i k i a [z,Sz}= I d 3 Xikl[z,Sz}, (4.83) 

where 5 Z describes a variation around the dyon satisfying the constraints. Indeed, in this 
case, k^ = = k e = k r , while k a = £ Qb . In other words k l = <pn is the electric 
potential on the horizon, while k 2 = —tpn is minus the magnetic potential on the horizon. 
Now, the results of the previous subsection imply that we get at infinity, 

/ d 3 x i k j u [z J 5z}= I d 3 x i k i u [z 1 z-z] = 5 z M- ( j) H 5 z Q-i/j H 5 z V. (4.84) 
For the matter part, we have 

d n ~ x x % C* M = ~ r dd r d<p e ab A a SB bLi , (4.85) 

S r 4?I " ^0 JO 
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which vanishes on the horizon r = r + where Aq vanishes. Note that in the gauge where 
Aq vanishes at infinity, the matter part gives no contribution at infinity, but 4>h8 z Q + 
i^hSz'P at the horizon, as it should. 

Finally at the horizon, the purely gravitational part gives 



d n - l Xl ki ra ^ = ^-5 Z A. (4.86) 

c 07T 



This can be shown for instance by using the fact that k^ rav ' 1 is the time-space compo- 

"d/dt 



nent of a conserved superpotential k^jl, that can be proved to coincide, for variations 



satisfying the linearized field equations and up to an irrelevant term of the form d^v^^, 
with the conserved superpotential considered in [|22l . In turn the latter has been shown to 
contribute -j^5 z A at the horizon. 

This concludes the geometric discussion of the first law 

S Z M = ^5 Z A + <j) H S z Q + ^ H 5 Z V, (4.87) 

for variations satisfying the linearized equations of motion around the Reissner-Nordstr0m 
dyon. 



4.8 Euclidean approach 

Our setup also allows us to complement the work of |]5||,ll6l by evaluating the partition 
function in the grand canonical ensemble, along the lines of [|23l . 

For the three commuting observables H, Q, P, we thus would like to compute 

Z[/3^ c ^ c ] = Tre"^-^-^ = e* G , (4.88) 

where ^g(P, ~ P4> C i ~P4' C ) iS m e Massieu potential for the grand canonical ensemble 
(see e.g. E4l in the present context), 

vM/3, -P<f>, -M) = S((H), (Q), <P» - /3(H) + (3<j>{Q) + M(P), (4.89) 
d^ G = -{H)df3 + (Q)d{p<t>) + (P)d(pij). (4.90) 

The path integral representation for this partition function is 

Z%<j>^\ = J V^e 1 * , (4.91) 

where $ represents all the fields (z A , u a ) together with appropriate ghost fields C a , C a 
E51 (see e.g. |(26| for a review). The appropriate action is, 

/ e T = J dr(i j d 3 x a A (z)d z A - (H - <\> C Q - ^ C P)) + ghost terms. (4.92) 
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The path integral is taken over all periodic paths in r with periodicity (3 and N — > 1, 

A — > 4> c , Z — > —ijj c for r — ► oo. 

We now notice that the transformation defined by 

7r ij -> — xtt^', Ai -> -zA, ^ -> -i^, iV* -> -iiV*, (4.93) 

with all other variables unchanged maps the action Ij to a real action when all (trans- 
formed) variables are real. The latter action differs from the Lorentzian action (14.821) by 
the fact that the terms involving and B tT in NH± and N l Hi have the opposite signs. 
For the purely gravitational part, this is as it should be in order that the path integral 
corresponds to one over Euclidean metrics after integration over the momenta 7r y '. 

The leading contribution to the path integral is given by the value of e 1 * evaluated at 
the classical solutions satisfying the specified boundary conditions, that is, (i) the fall-off 
conditions (|4.63I) - (I4.67I) . (ii) fixed values of the potentials (0 C , ip c ), and (iii) a fixed inverse 
temperature /3. The Reissner-Nordstr0m dyon (RND) described by the lapse and the 
spatial metric given in (|l.ll) and the matter fields (14.501) . (14.551) is such a solution if <p c = 
4>h an d ?p c = ipH with 4> H , ip H defined in (|1.6I) since all variables affected by the above 
transformation vanish in this case. Furthermore, this solution is time independent and 
satisfies the (modified) constraints so that I J reduces to surface integrals. For the matter 
part, we find directly that I™ at (RND) = (3(f)HQ+f3il>HP- For the gravitational part, it has 
been shown for instance in ll27l in the current Hamiltonian context that I§ rav (RND) = 
-f3M + \A, with A given in (11.51) . 

Assuming then that the dyon is the only extremum, it follows that to leading order, 

m G = —f3M +\A + /3(j> H Q + M H P, (4.94) 
which is the expected result. 

5 Conclusion 

In this paper we have generalized the manifestly duality invariant double potential for- 
malism [7], [8 J to include potentials for the longitudinal electric and magnetic fields thus 
turning the scalar potentials into non spurious Lagrange multipliers. By introducing ad- 
ditional pure gauge degrees of freedom on the classical level, which corresponds to an 
additional quartet [|28l on the quantum level, we have turned a topological conservation 
law, the magnetic charge, into a dynamical one. 

We have shown on the example of the Reissner-Nordstr0m dyon that the formalism 
is tailor-made for a treatment of black hole dyons by standard action based methods and 
allows one to compute in the grand canonical ensemble. How to explicitly resolve the 
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string-singularity of the Kerr-Newman dyon and derive its thermodynamics will be dis- 
cussed elsewhere. 

In our approach Dirac strings are only needed for the coupling to dynamical dyons and 
the derivation of the Lorentz force law. It would be interesting to understand whether there 
are applications of the formalism in the non-Abelian case or extensions to gravitational 
magnetic charge. 

Note: After the present work has been accepted for publication, references ll29l l30l 
have been called to our attention. In these references, a manifestly covariant double 
potential formalism is developed. As explained there, the problem is the occurence of a 
second "photon" that has to be removed in an ad hoc manner. 
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